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Abstract 

We study the decay of the metastable symmetric phase in the standard model 
at finite temperature. For the SU(2)-Higgs model the two wave function correction 
terms Z 9 (ip 2 ,T) and Z x (ip 2 ,T) of Higgs and Goldstone boson fields are calculated 
to one-loop order. We find that the derivative expansion of the effective action 
is reliable for Higgs masses smaller than the W-boson mass. We propose a new 
procedure to evaluate the decay rate by first integrating out the vector field and 
the components of the scalar fields with non-zero Matsubara frequencies. The 
static part of the scalar field is treated in the saddle point approximation. As 
a by-product we obtain a formula for the decay rate of a homogeneous unstable 
state. The course of the cosmological electroweak phase transition is evaluated 
numerically for different Higgs boson masses and non- vanishing magnetic mass of 
the gauge boson. For Higgs masses above ~ 60 GeV the latent heat can reheat the 
system to the critical temperature, which qualitatively changes the nature of the 
transition. 
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1 Introduction 



At high temperatures the spontaneously broken symmetry of electroweak interactions 
is restored This is a direct consequence of the Higgs mechanism of symmetry 

breaking, and the corresponding phase transition is an intriguing aspect of electroweak 
interactions. It is also of great cosmological importance since at temperatures of order 
the critical temperature of the phase transition baryon-number violating processes fall 
out of thermal equilibrium [|J. As a result, the present cosmological baryon asymmetry 
is finally determined at the electroweak phase transition. 

Theoretical descriptions of the phase transition start from the finite-temperature 
effective potential whose local minima yield the free energy of states with a homogeneous 
Higgs field given by the position of the local minimum (cf. Considerable effort has 



been devoted to evaluate the effective potential in perturbation theory |q|-pO[ where, 
due to infrared divergences, a useful expansion is only obtained after a resummation of 
an infinite number of terms. Alternatively, interesting first results have been obtained 
by lattice Monte Carlo simulations |TT|-Jl3[, reduced three-dimensional action [|TJJ] and 



by using an average action at finite temperature [15 



The decay of a metastable homogeneous state involves "critical droplets" and there- 
fore inhomogeneous field configurations. The free energy of these configurations is usually 
estimated based on an effective action which is the sum of the finite-temperature effec- 
tive potential and the canonical kinetic term of the Higgs field. We will improve this 
approximation by evaluating the finite-temperature wave function correction terms of 
the SU(2)-Higgs model to one-loop order. This will allow us to determine the range of 
Higgs boson masses for which the derivative expansion of the effective action is reliable. 

Following Langer's theory of metastability |1| we evaluate the decay rate in the 
semiclassical approximation where one starts from some approximate "coarse-grained" 
effective action and expands around its saddle point interpolating between the symmetric 
and the broken phase. For finite-temperature quantum field theories the "coarse-grained" 
effective action is frequently approximated by the one-loop effective potential together 
with canonical kinetic terms for the scalar fields. The Callan-Coleman theory for the 



decay of the false vacuum at zero temperature in four dimensions [IT] is then applied to 



the effective three-dimensional theory at finite temperature [18 



One problem, which points to some inconsistency of this approximation, is that the 
occuring effective potential is in general complex. The origin of this problem in the 
calculation of the decay rate is the incorrect treatment of the scalar fluctuations which 
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appear twice, as contribution to the effective potential and as ffuctuations around the 
saddle point. Instead, one has to split the functional integration into two parts, such 
that the first one yields a "coarse-grained" effective action which has a barrier between 
the symmetric and the broken phase. In the second part a perturbative expansion is 
then performed around the non-trivial saddle point, and only here long wave length 
scalar fluctuations appear. In this spirit, in JTPJ the integration over the vector fields has 
been performed in the first step. The same approach has been followed in f° r th e 
analogous problem of the Coleman- Weinberg model at zero temperature, in citebfhw for 
nonabelian gauge theories and in [21] for Yukawa models at finite temperature. Here 



we will further improve these calculations by integrating in the first step also over the 
components of the scalar fields with non-zero Matsubara frequencies. This will lead us 
to a "coarse-grained" effective action which is similar to, but not identical with the one 
used by Dine et al. |J. 

Several aspects of the cosmological electroweak phase transition, such as nucleation 
and growth of critical droplets, wall velocity, transition time etc. have already been 



discussed in the literature (cf. ||, ||,[ga]-[25|). We shall use our results on the finite- 



temperature effective action to study in more detail the dependence of the course of 
the phase transition on the Higgs boson mass and also the influence of a non-zero mag- 
netic mass of the W-boson. Our numerical calculations show that reheating effects can 
qualitatively change the nature of the transition for Higgs masses larger than ~ 60 GeV. 

The paper is organized as follows. In sect. ^] we evaluate the two wave function 
correction factors for the SU(2)-Higgs model to one- loop order. The decay of metastable 
and unstable homogeneous states is discussed in sect. ^. The course of the phase tran- 
sition is studied analytically and also in detail numerically in sect. f|[ The main results 
are summarized in sect. 5. 

2 The effective action at finite temperature 

Consider the SU(2)-Higgs model at finite temperature which is described by the action 
Sf,[$, W] = fdx (-Iw^W 9 ^ + (D^y(D^) - V (^ , (1) 



where 



dx = [ dr I d 3 x , (3 = i (2) 
p Jo J T 
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and W® u are the covariant derivative and the Yang-Mills field strength, respectively. 
The gauge fixing part of the lagrangian has been omitted, and we will always work in 
Landau gauge. The SU(2)-Higgs model contains the essential features of the standard 
model of electroweak interactions. Setting sinOw = 0, one neglects the mass difference 
between Z- and W-boson relative to the W-boson mass. This corresponds to the present 
accuracy in the description of the electroweak phase transition. We include the effects of 
the three generations of quarks and leptons with the usual gauge and Yukawa couplings. 

In order to find the thermal equilibrium states and to compute properties of the 
first-order phase transition we have to know the effective action rg[$] of the theory at 
finite temperature as functional of the Higgs field $ which, after gauge fixing, plays the 
role of an order parameter. The effective action can be expressed as a functional integral 

e -r/9[S] = jT[D6][D^ , (3) 

where the integration extends over fields with periodic boundary conditions in r. Again, 
gauge fixing and ghost terms have not been written explicitly. For stationary, r- 
independent fields the effective action corresponds to the free energy density of 

the system, 

T P [$\=PF[$,T\ . (4) 

Ordinary perturbation theory to any finite order does not yield a good approximation to 
r^f^] as the scalar masses are negative for small values of (p, and because of the infrared 
behaviour of the finite-temperature perturbation series. 

A systematic expansion in the gauge coupling g and the scalar self-coupling A is 
obtained by means of a resummation, where in all propagators the tree-level masses are 
replaced by one-loop plasma masses to order g 2 and A (cf. 0). The mass differences 
are then treated as counter terms, 

SSp = ~ J dx [W a » (8m\P Lia , + 6m 2 T P T , u ) W av - 5m 2 y - 5m 2 xX 2 ] . (5) 

Here Pl and Pt are the usual projection operators on longitudinal and transverse degrees 
of freedom of the vector field (cf. ||). The plasma mass terms 5m 2 are the differences 
between the masses 

m L = —g T + , (6) 

m T = —r~ > (<) 
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ml = (A^ + ^](T 2 -T 6 2 )+A^ , (9) 



and the tree-level masses 



m 2 = \g\ 2 , ml = \{?, V 2 -v 2 ) , m\ = \(<p 2 - v 2 ) . (10) 

In the plasma masses we have included the effect of the top-quark loop, with m t = 
ftv/\/2. The temperature at which the barrier between the symmetric and the 
broken phase vanishes, reads in leading order of the couplings 

2 lQXv 2 
Tb = 3g 2 + 8A + Iff ' ( } 

Note, that the tree-level Higgs- and Goldstone-boson masses are obtained from the scalar 
potential by (<f 2 = £f =1 fifi) 

^- rKPtj+KPh , (12) 

where 

are the projection operators on the Higgs and the Goldstone fields. 

At one-loop order the improved perturbation theory yields the effective potential to 
order g 3 , A 3/2 , 

Veff^ 2 ,T) = \{^ + \ + \fl)(T 2 -T 2 ) V 2 + \^ 



-(3ml + 6m* + m% + Sm 3 )^ + 0(<j s . A 2 . //). (11) 



which is equivalent to the result of the ring summation |7j . Several two- loop contributions 
are also known. For the SU(2)-Higgs model the effective potential has been calculated 
to order g 4 , A B, and for the abelian Higgs model it is known to order e 4 , A 2 |26l . 



The strength of the electroweak phase transition is rather sensitive to the magnetic 
mass of the gauge bosons, whose calculation requires non-perturbative techniques. In 
Landau gauge the one-loop gap equations yield m T = (g 2 /3n)T at <p = P|,PD|. In 
order to estimate its effect on parameters of the phase transition we will replace eq. (|7|) 
by 

2 2 9 rr2 i 9 2 / 1t -\ 

m T = l 7T^ T + ~rV , 15 
97T 4 



5 



and compute sensitive quantities for different values of 7. 

The effective potential flT3|) has degenerate local minima at tp — and tp = <p c > 
at a critical temperature T c and therefore predicts a first-order phase transition. The 
evaluation of the transition rate requires knowledge of a stationary point of the free 
energy which interpolates between the two local minima. The effective action (^) can be 
expanded in powers of derivatives, and for time-independent fields one has 

r>[$] = (3 J d 3 x (V e //(^ 2 , T) + hfu + Zjj($, T))V(fiV(fj + . . .) . (16) 

From eqs. @ and (|16"D it is obvious that the functions Zjj can be obtained from the 
inverse scalar propagator in the homogeneous scalar background field $ with spacial 
momentum k, 

dtpidtpj 

+(5 I j + Z IJ ($,T))k 2 + 0(i A ) . (17) 

The one-loop contributions to Zjj are shown in fig. (TI]). It is a well-known problem in 
finite temperature field theory that the order of a given Feynman graph can be smaller 
than the order coming from the vertices. This is due to the fact that loop diagrams ac- 
cumulate couplings in the denominators which decrease the power of the usual couplings 
coming from the Feynman rules. In appendix A we show that only graphs of the type 
shown in fig. 1 give contributions to the wave function correction in leading order, where 
one has to include not only the tree level but also the plasma masses for internal lines. 
Note, that the knowledge of the effective action for static fields is not sufficient if one is 



interested in time-dependent processes such as the thermalization of perturbations |27 

Eq. (|16|) is an expansion in powers of (V/m^) 2 , where denotes all masses which 
appear in propagators of internal lines. Clearly, the terms with zero Matsubara frequen- 
cies for the internal lines dominate, since they have the smallest denominators. For the 
same reason the top-quark contribution is suppressed, since all Matsubara frequencies 
are non-zero for fermions. It is straightforward to evaluate the contributions of fig. ([!]) 
as function of the couplings g and A, the scalar propagators 

A(k 2 } m*) = 7 -^— , (18) 
k z + vrq 

and the vector propagator 

D^(k) = /\(k\ml)P L ^ + /\(k\m 2 T )P T ^ . (19) 
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After some algebra one obtains a sum of integrals of the type 

d 3 l . 



T 



T 



47T | fc 
1 



(2tt) 3 
arctan 



A(/ 2 ,m 2 ) A((/ + £;) 2 ,m 2 ) 

ife i 



T 



mi + m2 
( 1 



A- 2 



47rmi+m2 \ 3 (mi +m2) 2 
Collecting all terms we obtain the final result 

Zij($, T) = Z^ 3 , T) Jft + Z x (<p 2 , T)Pfj 



+ 0( 



(20) 
(21) 

(22) 
(23) 



where 



2 Am 2 I _ + _ ) - 2g< 



mr, 



i 2 2 

+ -<? m 



1 5 



m x + m T 



2m| 



m 



T , 



and 



2Am 2 



1 
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(24) 



(25) 



6n [(m^ + m x ) 3 " \m x + rriT vn 9 + tut , 
with m 2 = \ip 2 and m 2 = g 2 (p 2 /4. In terms of the complex field $ this result can be 
written as 

^jVyp/V^j = ^($ t $)- 1 [V($ t $)] 2 ^+fw t V$ - ^($ t $)- 1 [V($ t $)] 2 ) Z x . (26) 



(For completeness we give in appendix A the wave function correction term Z((p 2 , T) 
for the abelian Higgs model.) The Z-factor for the Higgs field is dominated by the 
contribution from the transverse vector boson loop which, if one neglects the magnetic 
mass, diverges at if ~ 0, 

5g 2 T 



' 16n tp 

Neglecting plasma mass terms for the vector bosons and using tul = 
obtains instead 

y f 2 ^ 21^T 

Z ^' T) ~"32^ ■ 

This agrees with the result of a previous calculation by Moss et al. 

different technique. The comparison of eqs. ( |27| ) and (^) shows the important difference 

between the ordinary and the improved perturbation theory with plasma massesQ The 

: In a very recent paper |29 ] Z v has also been calculated using the improved perturbation theory. Our 
results differ from those obtained by these authors. The differences can be traced back to the treatment 
of the mass counter terms of eq. ([s]) . 



(27) 

rriT 3> m x one 
(28) 

who used a 
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order of magnitude for Z^ is the same in both cases, the sign, however, has changed. 
Wave function correction terms in the 't Hooft-Feynman background gauge have recently 



been evaluated in |MJ. Note, that despite the divergence of Z^ at ip ~ the correction 



to the surface tension o = Jq c dip y 2(1 + Z^ip 2 , T c ))V e ff(ip 2 , T c ) is finite. 

What is the effect of the wave function factors Z(ip 2 ,T) on parameters describing 
the electroweak phase transition? The first-order phase transition is due to the barrier 
in the effective potential ( JH| ) which is essentially generated by the vector-boson loops. 
This implies that at the critical temperature T c the magnitude of the local minimum in 
the broken phase is ip c /T c ~ g 3 /X- One also easily verifies that (T 2 — T 2 )/T 2 ~ fi ,4 /A. 
Hence, in the effective potential fll4[) the terms quadratic, cubic and quartic in ip are 
all of order ~ (g 12 /A 3 )T c 4 . The perturbative approach is only reliable if this energy 
density is larger than the density due to the magnetic mass of the W-boson, which is of 
order m\T c ~ g 6 T 4 . Hence a necessary condition for the smallness of non-perturbative 
effects is A < g 2 . Further, if one uses the high-temperature expansion one has to satisfy 
gip c /T c < 1, i.e., A > g 4 '. In this case the allowed range of A is given by g A < A < g 2 . 

A naive estimate of the size of Z(f 2 ,T c ) can now be obtained based on eqs. @ - 
(|9]), ( p4[) and (p5|) . One easily verifies that, at ip c and T c , the plasma masses are of order 

4 3 

m T „ —T c < m L , m VtX ~ ~^ T c ■ ( 29 ) 
From eqs. (|2^), (|25| ) one then concludes 

Z<px = °(tz) ■ (3°) 



Hence, for the allowed range of couplings determined from the effective potential correc- 
tions due to the Z-factors can be treated as perturbation. 

Clearly, these estimates are rather naive. As the saddle point of the effective action, 
which interpolates between the two local minima, varies from ip = to ip = p c , the full 
scalar mass terms 



2 _ d 2 V eff (ip\T) 2 _l dV eff (ip\T) 
m *~ df~ 2 ' m *~f df~ ' 



become negative and the expansion in powers of derivatives breaks down. Even for 
positive mass terms this expansion is problematic for stationary points of the effective 
action. For these field configurations 

VV = ^V eff (ip\T) 

= m 2 x ip. (32) 
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the equation of motion shows that V/m* = 0(1). Hence, higher order terms in the 
derivative expansion are non-negligible As we shall see in the following section, a more 
careful treatment nevertheless confirms the above conclusion on the size of the wave 
function correction terms. 



3 Decay of the false vacuum 

At temperatures below the critical temperature T c the symmetric phase becomes 
metastable and decays via nucleation and growth of droplets of the broken phase. In 
condensed matter physics a theory for the decay of such metastable states has been de- 



veloped by Langer many years ago |fLq| . Starting from the Fokker-Planck equation for 
the probability distribution of droplets with different sizes, Langer derived an expression 
for the nucleation rate in terms of a saddle-point (p of the free energy which interpolates 
between the two phases. A similar result for the decay rate has been obtained by ap- 
plying the Callan- Coleman theory for the decay of the false vacuum at zero temperature 



in four dimensions [|17|] to the effective three-dimensional theory at finite temperature 
[ p~8] . We will base our discussion on Langer's theory according to which the decay rate 
is given by 

where Y p is given by eq. The "dynamical factor" k has recently been evaluated in 
terms of viscosities of the electroweak plasma f25 |. The rate (|33|) is determined using 



a time independent solution of the Fokker-Planck equation. This is justified because 
the characteristic time scales of the microscopic processes (e.g., WW scattering) are 
several orders of magnitude smaller than the time needed for the phase transition. Thus, 
the Langer theory is applicable, and it is sufficient to use a stationary solution of the 
Fokker-Planck equation. 

The non-trivial phase structure of the Higgs models at finite temperature becomes 
only apparent if quantum fluctuations are taken into account. The tree-level potential 
possesses only a single local minimum at (p ^ 0. This is analogous to the case of Coleman- 
Weinberg symmetry breaking by radiative corrections. In order to obtain the decay rate 
one first has to integrate out degrees of freedom which generate the barrier between 
the symmetric and the broken phase and only then, in a second step, determine the 
interpolating saddle point. In |19| it has been suggested to first integrate out the vector 



field. The same route has been followed in the abelian Higgs model at zero temperature 
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and Coleman- Weinberg symmetry breaking [20]. This method has then been extended 
to the standard model and Yukawa models [^j] . 



Here we suggest to improve this method in the following way. In eq. (|33|) the back- 
ground field $ depends only on spatial coordinates. Therefore in the classical action Sp, 
which appears in the integrand of the functional integral, only the quantum field $o is 
shifted, i.e., the mode with zero Matsubara frequency, 

$(r,f) =J2K(x)e tuJnT , u n = 2nnT . (34) 

n 

This suggests to integrate out the scalar fields with non-zero Matsubara frequencies in 
addition to the vector fields. Since the " Compton-wave lengths" of these modes are 
smaller than the wall thickness of the critical droplet, this procedure should yield an 
appropriate effective action which plays the role of the "coarse-grained" free energy in 
Langer's calculation. The scalar part of the functional integral is therefore split, 

D$ = ]\D$ n = D$ l[D$ n , (35) 

n n ^0 

and in the first step all modes except $o are integrated out. 

The functional integral over vector fields and the non-zero frequency modes of the 
scalar field yields an effective action Tp which replaces the classical action Sp in eq. (|33|), 

e~ T ^ = J^[D%] n / [D& n ] JjDW\ exp {-Sp[Q + $,W] + J d^^|U(x)) 

= ^[D$ ] exp ^-f>[$ + $ Q ] + J rf 3 x^^$o(x)^ . (36) 

At one-loop order one obtains for the effective action Tp a sum of terms of the type 
calculated in sect. |2|, where now the zero- frequency contribution of the scalar fields are 
absent in the loop integrals. In addition there is a non-local contribution which stems 
from the integration over the vector field, 

r„[*] = f«[$] + f( 2) [$] , 

tf[$\ = pjd 3 x(v eff (ip^T) + ^(5u + Z I j^ 2 ,T))V^ I V^j + ..}j , 
rf[<$>] = ^-(3 J d 3 xd 3 x'J^(x)D^(x,x')J^(x') , (37) 



where = ^(&d ll T a $ — 9 M $V a $) is the SU(2)-current and D T ^ V is the propagator for 
the transverse part of the vector field in the background field $(x). 
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The difference between the effective potentials is given by 



9^V 6// e// V ^ 7 (27r) 3 \fc2 + m 2 jfe2 +m 2^ 

3A 

= — p>T{m v + m x + C) . (38) 

Here C is a linearly divergent constant which vanishes in dimensional regularization. 
We have explicitly checked that for other regularizations it cancels against the sum 
of divergent terms which arise from integration over fields with non-zero Matsubara 
frequencies. To leading order in the couplings the masses m v and m x are given by eqs. 
(|D and (|). Integrating eq. fl38|) gives precisely the cubic scalar terms in the potential 
(fn|). Hence, defining the split ( |34"D of the functional integral in dimensional regularization 
one obtains to order 0(} 3 , \ 3 ^ e ) 

Veff( V 2 ,T) = 



where 




-(3ml + Qm 3 T )-?- + 0(g\\ 2 )T b * , (39) 

1Z7T 



11 2^2 , 9 2{ P 2 



ml = =-g' z T z + 



6 ° 4 

2 2 3 rr2 , 9 2 

The effective potential V e ff is similar to the potential used by Dine et al. H] to describe 
the electroweak phase transition. For (p <C T the term linear in T essentially becomes 
the one-loop vector boson contribution reduced by a factor 1/3 0. To estimate the 
size of non-perturbative corrections we have kept a non-zero magnetic mass with an 



order of magnitude given by the one- loop gap equations @,|T0|- A novel feature of the 
potential (|39|) is the way in which contributions from scalar loops are treated. There are 
no terms cubic in the scalar masses, since the integration over $ has not been carried 
out. However, contrary to the treatment in 0, the scalar-loop contributions to the term 
quadratic in (p are kept, since they arise from integration over $ n , n ^ 0. This has an 
important effect on the difference between critical temperature and barrier temperature 
of the transition. For m# = 60 GeV we find an increase by a factor of two. 

We conclude that keeping the contribution of scalar loops to the T 2 -term in the 
effective potential and dropping the cubic scalar mass terms yields the correct effective 
potential for the zero-Matsubara frequency part of the scalar field. 
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As discussed in sect. [|, the part of the Z-factors arises from the zero- frequency 
contributions to the loop integrals. Hence, pure vector boson loops give the dominant 
contribution to Z, since the integration over $o is not carried out. The result is easily 
read off from eqs. (|24j),(p5[), 

ZJip 2 ,T) = — g 2 m 2(J_ + JL) 
^ ' ; 32vr \2m 3 L m\) ' 

Z x (ip 2 ,T) = . (40) 

The "coarse-grained" effective action Tp is now given by eqs. (|3T|) , fl39|) and ([4(]). Note, 
is positive. Hence, the surface tension and therefore the strength of the phase transition 
is increased. Z^ is shown in fig. (0) for three different values of the magnetic mass. 

Using eqs. (|33|) - (|40"D we can now calculate the decay rate of the metastable sym- 
metric phase. As discussed in sect. 2, the Z-i actors Z^^ are expected to yield only small 
corrections. Therefore we start from Z = and the effective potential (|39D and determine 
a saddle point (p which interpolates between the symmetric and the broken phase. The 
remaining $ -i n t e g ra ti° n is then carried out in the Gaussian approximation. The spec- 
trum of fluctuations contains one negative eigenvalue A_ which guarantees that Zp[(p] is 
purely imaginary. Further, there are six zero modes from translational invariance and 
global SU(2)xU(l) symmetry, and a discrete and continuous positive spectrum. The 
contribution of the zero-modes to the fluctuation determinant is given by || 

/ 1 ~ - \ 3/2 

V trans V = V, 

V mt = ^(^W) 3/2 , (41) 

where V is the total volume of the physical three-dimensional space. In the thin wall 
approximation also the contribution from "Goldstone modes" has been evaluated 0, 
which correspond to deformations of the droplet surface. However, for the electroweak 
phase transition the thin wall approximation is known not to be reliable @,[25|. We 
therefore base our calculations on the decay rate (cf. (|33|) , (|4l|) ) 

£ = K | A_ |"Va v uV h\/ V)\ u ^}-^ ] = A r®P$] , (42) 

where (p is the saddle point and ji = m^OjT). The negative eigenvalue A_ is approxi- 
mated by — 2/R 2 , its thin wall value, where the radius R of the critical droplet is defined 
by d 2 ip(r = R,T)/dr 2 = 0. k can be estimated using the results given in ||25|| . 

We have numerically determined the saddle point (p at temperatures close to the crit- 
ical temperature T c for different values of the scalar self-coupling A and the parameter 7, 
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i.e., the magnetic mass mx- This will be described in more detail in the following section. 
Fig. (^|) shows the ratio of the pre-factor in eq. fl42|) normalized to T 4 as function of the 



Higgs mass at the respective temperature T e where the phase transition is completed. 
Clearly, the logarithm of the pre-factor is much smaller than Tp[0\ which, at the tem- 
perature T e , is ~ 140. Hence, the semiclassical approximation is justified. The difference 
between the pre-factor A and T 4 is numerically unimportant for the description of the 
cosmological phase transition. As the Higgs mass m# approaches the critical mass m c ^ lt , 
where the transition becomes second-order for 7 ^ 0, the pre-factor strongly decreases. 
This is mostly due to the decrease of fi. 

A measure for the size of the one-loop correction to the Z-factor is the ratio 



(43) 



This correction is shown in fig. as function of the Higgs mass for different values of 
7. For 7 = the perturbative expansion breaks down at tuh ~ 80 GeV, in agreement 
with results obtained previously ||. For 7 = 1, 2 the correction never exceeds 10% which 
demonstrates the importance of the infrared cutoff provided by the magnetic mass. 

The effective potential to order 0(} 3 , A 9 / e ) is given by eq. ([14]). The corresponding 
quantity without the one-loop contribution of the static part of the scalar field <3? is 
V e ff (|39|). A more accurate approximation to the effective potential can be obtained by 
performing the ^-integration over the expression (^) in the Gaussian approximation. 
The result is the potential 

Veff = V eff -^{m% + 3m 3 x ) , (44) 
where, contrary to eq. (|14"D, the scalar masses are now given by 

ml=^ff\ ml = . (45) 



These masses agree with the mass terms (§) and (|9|) only to leading order in the couplings. 
The potentials V e ff and V e ff clearly differ by an infinite number of terms, if the masses 
(f45|) are formally expanded in powers of A. The most important property of the mass 
terms (^) is that they become negative as the scalar field varies between the two local 
minima of V e ff- Hence, the true effective potential V is complex. This general feature of 
non-convex potentials, which is well known for models at zero temperature, is not visible 
to any finite order in perturbation theory. 
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What is the meaning of the imaginary part of the effective potential? At zero tem- 
perature the imaginary part is defined by first substituting 

m 2 (tp) — > m 2 ({p) — ie , (46) 

and then performing the limit e — > 0. This yields the imaginary part 

I(<p) = 2 hnV(<p) (47) 
= I rn 2 (ip) | 2 Q(-m 2 ((f)) (48) 

which is interpreted as decay rate per unit time and volume of the homogeneous state 
into inhomogeneous states of lower energy |31j . 

It is straightforward to adopt the same procedure at finite temperature. One imme- 
diately obtains from the cubic scalar mass terms in eq. (f4"4"D 

Tfa) = 2 ImVefffaT) (49) 

= ^(l<l 3/2 eK)+mi 3/2 eK)) ■ (so) 

This should then be the finite-temperature decay width of a homogeneous state into 
mixed states. This interpretation appears sensible at a stationary point (p of the potential 
V(f 2 ,T) where m x = 0, i.e., eq. (|4"9| ) should be applicable at the top of the barrier. For 
arbitrary values of cp the interpretation of the rate ( f4"9|) is less clear due to the gauge 
dependence of this expression. 



4 The cosmological phase transition 

We can now use the results obtained in the last section to investigate the cosmological 
electroweak phase transition. This phase transition has already been studied in consid- 
erable detail by several groups P1,P|,P2"1-|2~5|. Here we will extend the previous studies 
mainly in two respects. We shall study the dependence of the transition on the mass of 
the Higgs boson, and we shall also investigate the influence of a non-zero magnetic mass. 
In this case the transition becomes second order for Higgs masses above a critical mass 
m c ^ lt 0, and the reheating of the symmetric phase can qualitatively change the nature 
of the transition. 

The main features of the transition can be understood by considering the simplified 
potential 

V= l -a{T 2 -T 2 ) V 2 - l -bT^ + \\^ , (51) 
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where the parameters a and b can be read off from eq. (0). Here we neglect tul, which 
contributes essentially only to the T 4 -term, and set 7 = 0. An analytical treatment 
of the transition is possible in the thin wall approximation Quantitatively, this 

approximation is too crude for the electroweak phase transition. However, it is still 
useful to get a qualitative picture of the transition and to obtain the order of magnitude 
of the various parameters involved. This will be confirmed by our numerical calculations 
discussed below. 

In the thin wall approximation the nucleation rate is given by 

m_ Aew (_ B( t^y\ _ (52) 



V \ 1 t, - f 

where 



A = uTb ' B = 3^o a 2 W 2 ' (53) 
with uj = 0(oo). Tf, is the barrier temperature, and t c is the time at which the critical 
temperature T c « T b is reached. Using eqs. (0), (|39|) and the relation between time and 
temperature, t « 0.03 mpi/T 2 with mpL = 1.2 • 10 19 GeV [3^], one obtains the relations 

T b = 120GeVG F 1/2 (3m^ + m^ + 2m^)- 1 /^_^_) , (54) 

t c = 1.7 ■ 10~ n sec G F (3m 2 w + m 2 H + 2m 2 ,) ^ 100GeV N ) (55) 

V m H J 

9 10 15 

D ~ 04 _4 /0__2 , __2 , r> 2N9 .__7 • 



3 4 vr 4 (3m^ + 771^ + 2m t 2 ) 2 m 7 H 



(57) 



From these equations one can read off the involved orders of magnitude of time and 
temperature and also the dependence on the Higgs boson mass. 

As the universe expands the temperature T decreases. Once the critical temperature 
T c is reached critical droplets can start to nucleate and grow. At time t the number 
density of droplets of size r is 

n(t,r) = vU-^j . (58) 

Here v is the velocity of the droplet wall, and the initial size of the critical droplet has 
been neglected. The fraction of space filled with the new phase is then 

/ 4-n- rr ma x(t) \ 

f B (t) = exp [ -— ^ drr 6 n(t,r)\ , r max (t) = v(t - t c ) , (59) 
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where the exponentiation in eq. (|59"D accounts for the overlap between different droplets. 

The phase transition is completed at time t e , which one may define implicitly by 
/s(^e) = 1/e. The nucleation rate is non-negligible only for times very close to t e . This 
allows one to derive simple expressions for the time t e , the "transition time" r during 
which the phase transition essentially takes place, the total number of droplets N(t) per 
unit volume and the average radius R(t), which are defined as 

rr m ax(t) 

N(t) = / drn(t } r) , (60) 



o 



R(t) = — |— r max(t) dr r n(t,r) . (61) 



N(t) Jo 

A straightforward calculation yields for the time t e the implicit relation 



and for the transition time r one obtains 



IBtl 



(63) 



The fraction of converted space, the number of droplets and the average radius are given 
by 

f B (t) « expf-expf^)) , (64) 



T 

N(t) « J_exp(^] . (60) 



■57TT 3 V T 



R(t) w r . (66) 

These equations demonstrate that the phase transition takes place in the time interval 
t before t e . Using eqs. (|55| ) , (|56|) , (|62|) and (|63D one easily verifies the following hierarchy 
between the different time scales 

t c > t e - t c > t . (67) 



We have also studied the phase transition numerically. Here the free energy of the 
critical droplet is calculated using the potential V e ff (cf. ([39])) derived in sect. |[ We keep 
a non-zero magnetic mass of the size given by the one-loop gap equations, i.e., we choose 
7 = 1 in eq. (|T3|). For the potential (^) the stationary point (p and the corresponding 
free energy are determined numerically for temperatures below T c down to T e , where 
the transition is completed. The course of the transition depends on the velocity of the 
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droplet wall which has been discussed in detail in the literature P|,f2"3fl,[pl|]. Here we just 
illustrate the dependence by calculating the course of the transition for three different 
velocities, v — c, v — c/10 and v = c/1000. The initial radius of the critical droplet, 
which is very small (cf. ||), is also taken into account. 

Figs. ([5p - (|5p show the result for m# = 30 GeV, which we have chosen as an example 
of a small Higgs mass. The completion time t e is ~ 3 ■ 10 -12 sec and depends only weakly 
on the wall velocity (fig. (|)). As one may have expected, the number of droplets (fig. 
(|5|)) and the average radius (fig. (|^)) at the end of the transition scale approximately 
like v~ 3 and v" 1 . The number density of droplets is given per cm 3 . Note, that the size 



of the horizon at the critical temperature is of order 1 cm. Following Enqvist et al. [22 



and Carrington and Kapusta we have plotted in fig. (§) the change in temperature 
as function of time. The system supercools from time t c until the time t s at which the 
supercooling stops and the liberated latent heat leads to some reheating. The increase 
in temperature can be estimated by 

AT R = - L , (68) 

c v 

where the latent heat is given by the derivative of the effective potential, 

d ~ 



(69) 

T=T C 



and the specific heat of the symmetric phase is 

Cv = jL e (T) , e(T) = — vr 2 T 4 . (70) 
dT K 1 ' v ; 120 v ; 

Here e(T) is the energy density of the standard model at high temperatures to leading 
order in the coupling constants, where gluons and the U(l)-gauge boson have been 
included. The above equations are easily derived from the usual thermodynamic relations 
and the relations for the pressure p s (T) = —V e ff(T,0) and Pb(T) = —V e ff(T,ip 2 (T)) in 
the symmetric and broken phase, respectively. Using these equations to compute the 
reheating temperature assumes that the wall velocity is sufficiently small so that the 
latent heat can indeed be converted into thermal energy of the symmetric phase. 

In fig. (|) the true evolution is compared with a transition where thermal equilibrium 
is maintained at all times. In this hypothetical case the temperature remains at T c until 
the time t' e , at which all latent heat has been used up to expand the universe. At later 
times the cooling due to the expansion then proceeds in the usual manner. Note that 
the real phase transition starts after considerable supercooling and completes in a very 
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short time. The temperature after reheating is slightly larger than the temperature of 
the adiabatic transition at the same time. This is due to the entropy generated in the 
transition. 

For the Higgs mass m# = 60 GeV the transition looks qualitatively the same. The 
time interval t s — t c decreases by about two orders of magnitude to ~ 2 • 10~ 14 sec (fig. 
(|9|)). Correspondingly, the average radius decreases by about two orders of magnitude 
whereas the number density of droplets increases by about six orders. Since much less 
time is available for supercooling, the temperature after reheating is now closer to the 
critical temperature T c than for m# = 30 GeV. Our numerical results for m# = 60 GeV 



are consistent with those obtained in ||25|| , given the differences in the choice of m t , the 
magnetic mass and the wall velocity. 



A dramatic change occurs for m# = 80 GeV (fig. (10)). After the initial supercooling 
period the liberated latent heat starts to reheat the system at the time t s < t' e , long before 
the transition is completed. Due to the small supercooling the latent heat can reheat 
the system to the critical temperature. After that no new droplets nucleate, the average 
radius R increases linear with time, and the transition completes in equilibrium (fig. 
fllDD), resembling the QCD phase transition [22|. Compared to m# = 60 GeV the time 
interval t s —t c has decreased again by two orders of magnitude to ~ 3- 10~ 16 sec. Contrary 
to the QCD phase transition the time difference t' e — t c is much smaller than the Hubble 
time due to the smallness of the latent heat. 

The change of the nature of the phase transition for large values of the Higgs mass 
is a consequence of the magnetic mass. For any finite value of this mass, the phase 
transition becomes very weakly first-order with rising Higgs mass and eventually second- 
order above a critical Higgs mass m% %t || . Our calculations show that below this critical 
value there is a large range of Higgs masses for which the universe reheats up to T c . 
This is shown in fig. (|Tl|) where the reheating temperature ATr is compared with the 
supercooling temperature needed to initiate the phase transition AT S c- Clearly, the 
reheating temperature reaches T c if ATr > AT$c- 

This interesting result can also be understood analytically. Consider a potential with 
a generic plasma mass, in our case the simplified potential fl5T|) with non-zero magnetic 
mass, 

V = ^a{T 2 -T^ 2 - l -bT{c 2 T 2 + V 2 fl 2 + -^ . (71) 

This potential implies a change from a first-order to a second-order phase transition as 
the Higgs mass increases. Since there are strong hints that the magnetic mass is non- 
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zero, we expect that eq. ([71]) is a reasonable approximation for the effective potential of 
the standard model. Above a critical value of the scalar coupling, 



A 



crit 



b 

2c 



(72) 



the potential Q7Ip describes a second-order phase transition. For values of A slightly 
below X cr it the transition is weakly first-order. We are interested in the behaviour of 
several quantities as A approaches X cr u from below. The position ip c of the degener- 
ate minimum at T c vanishes in this limit. Some other characteristic quantities of the 
first-order transition vanish as powers of ip c . We find for the mass difference Am, the 
maximum of the effective potential V max , the surface tension a, the transition rate T, 
the latent heat L and the temperature differences T c — T& and T c — T s , where T s is the 
temperature corresponding to t s , 



Am 



m 



Crit im 

H ~ m H 






ymax 


r*j 




a 


r\j 




L 


rsj 




r 






T c -T b 


rvj 




T c -T s 







(73) 



Using these equations one obtains for t s , the time at which the supercooling stops for 
the true transition, and for t' e , the time for which the equilibrium transitions complete, 
respectively, 



t - tr 



L 



4 

Am 



Am' 



(74) 



Clearly, for Higgs masses very close to m c ff l one has t s < t' e , i.e., the supercooling 
ends, the liberation of the latent heat reheats the system to T c and the phase transition 
completes in equilibrium. For vanishing magnetic mass this phenomenon does not occur 
for Higgs masses smaller than the vector boson mass, for which the perturbative approach 
is reliable. 
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5 Summary 



In this paper we have extended previous studies of the electroweak phase transition in 
several respects. In particular, we have studied in greater detail the decay rate for the 
metastable symmetric phase and also the course of the phase transition. On the whole 
the by now familiar picture is confirmed that the phase transition is weakly first-order for 
Higgs boson masses up to the W-boson mass where the perturbative approach becomes 
unreliable. 

We have calculated the two wave function correction factors Z v and Z x for the SU(2)- 
Higgs model, and we have evaluated the corresponding correction to the free energy of 
critical droplets. We find that the perturbative expansion is reliable for Higgs masses 
below 80 GeV. 

The evaluation of the decay rate of the metastable symmetric phase is a non-trivial 
problem, especially since the metastability is entirely due to quantum corrections. Hence, 
two expansions are needed, first an improved perturbative expansion to obtain the 
"coarse-grained" effective action with a barrier between the two local minima, and sec- 
ond the expansion around the critical droplet, the saddle point of the "coarse-grained" 
effective action. We have improved previous calculations by integrating out the vector 
fields and the components of the scalar fields with non-zero Matsubara frequencies in 
the first step. However, further work is still needed to obtain a better understanding of 
renormalization, gauge dependence and convergence of the performed expansions. 

As a by-product of our calculation we have derived a formula for the decay rate of 
a homogeneous unstable state, which is the finite-temperature analogue of the familiar 
zero-temperature decay rate. Both formulae are obtained from the imaginary part of the 
effective potential. The interpretation of the formula away from a local maximum of the 
potential, in particular its gauge dependence, requires further study. 

Our numerical analysis of the course of the cosmological electroweak phase transition 
shows that at a Higgs boson mass of ~ 60 GeV a qualitative change of the nature of 
the transition occurs. This is a consequence of the very weak first-order phase transition 
which in our case is due to the assumed non-zero magnetic mass of the W-boson. For 
Higgs masses above 60 GeV the reheating temperature reaches the critical temperature 
before the transition is completed. The phase transition then continues in equilibrium, 
similar to the QCD phase transition. This may have interesting implications for the 
generated density fluctuations. 
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As the Higgs mass increases the phase transition becomes more and more dependent 
on non-perturbative properties of the symmetric phase. In our calculations this is re- 
flected in the strong dependence on the value of the magnetic mass of the vector boson. 
Recent lattice calculations indicate that such non-perturbative effects may increase the 
strength of the first-order transition. A deeper understanding of these non-perturbative 
effects appears to be a crucial step on the way towards a theory of the electroweak phase 
transition. 

The work of D. B. has been supported by the "Graduiertenkolleg fur theoretische 
Elementarteilchenphysik" , Universitat Hamburg. Z. F. acknowledges partial support 
from Hung. Sci. Grant under Contract No. OTKA-F1041. 
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Appendix A Contributions to Z($,T)) 

As it was promised we show that only graphs of the type shown in fig. 1 with plasma 
mass insertions on the internal lines give contributions in leading order to Z($,T). 

Instead of g, V\ and f t we use in this appendix a generic coupling h. 

A diagram can contribute to an order in h lower than the power in h obtained 
from vertices, if some of the Matsubara frequencies of the loop variables vanish. In the 
following we call these loops and the corresponding lines soft . Consider a contribution 
to the effective potential from a graph G of the following type: 

a. it has at least one internal scalar line I scalar'-, 

b. by cutting l sca i ar the diagram remains one-particle irreducible; 

c. r loop variables (among them the momentum of l SC aiar) are soft; 

d. the loops of non-zero Matsubara frequencies (in the following hard loops) yield leading 
order plasma mass corrections (cf. fig. 12, where thick lines are hard, thin ones are soft, 
r = 2). 

According to appendix A of this graph gives a contribution Vg to the effective 
potential of order h r+2 . Writing Vg as a Feynman integral over p (the momentum of 
lscaiar) with the appropriate self-energy insertion S G (p = O,^ 2 ) and performing a Taylor- 
expansion in p 2 one gets 

Vg = I it P \3 -2 Z — 2 j: G(po = 0,f) 
J (zir)-* p z + m z 

= I ft\ -2 1 2 \^g(Po = o,f = o)+fz G ( Po = oy = o)+ ..." 

J {ztt) 6 p z + m z 1 J 

= C + h3 J^^l^ Po = °^ = 0) > (A ' 1} 
where H' G = dUc/dp*, and we have used the substitution p = mx (remember that m 
is of order h). Since Vq is of order h r+2 and the above integral over x is independent 
of h, the derivative term Ti'(p = O,^ 2 = 0) is of order h r ~ l . In our case one-soft-loop 
contributions with plasma mass corrections to S(po = O,^ 2 ) correspond to two-soft- 
loop (r = 2) contributions to Vq, thus the wave function correction term is at least of 
order h 2 ~ l = h. Summing all plasma mass corrections yields the full propagators of the 
improved perturbation theory for internal lines. This 0({) contribution is non-zero and 
it is given by eqs. (Pi|) , (P5|) . 

According to (A. 16) of other diagrams (eg., diagrams with more soft loops or more 
complicated hard loop structures then mentioned above) lead to higher order corrections 
to Z($,T). 
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For completeness we also list the Z-factors for a model with n complex fields 
$ = ($ (1) ,...,$ (n) ), = (v?! a) +iy4 a) )v / 2~, where the diagonal global U(l) symme- 
try is gauged. At one loop order, zero-Matsubara frequency internal lines yield the wave 
function correction terms (cf. fig. 1. a,b): 

lzlf>(*,T)V<p<T ) V<pV>= (A.2) 
- i($ t $)~ 1 (V$ t $ - $ f V$) 2 Zi + V$ f V$Z 2 + i($ t $)~ 1 (V($ t $)) 2 Z3, 



where 



Z 1 (<p 2 ,T) 



2e 2 T 



3/T \ to v + tot + m ^ / 

2e 2 T 1 



37T TO X + TO j 1 



„ , o _ x e 2 Tm 2 / 1 10 \ /a \ 

e is the gauge coupling and to = &p. For n = 1, the usual abelian Higgs model, only two 
linear combinations contribute, Z x = Z 2 — Zi and Z v = Z 2 + Z 3 . 
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Figure captions 

Figure [I]: One-loop contributions to the wave function correction factors. 

Figure The wave function correction factor for three different values of the mag- 
netic mass. 

Figure |3]: The pre-factor of the decay rate as function of the Higgs boson mass. 

Figure ^]: One-loop wave function correction as function of the Higgs boson mass itlh 
for different values of 7. 

Figure Volume fraction of the new phase as function of t — t c ; itlh = 30 GeV. 

Figure |6]: Number density of droplets per cm 3 as function of t — t c ; m H = 30 GeV. 

Figure |7|: Average droplet radius as function of t — t c ; itlh = 30 GeV. 

Figure |8]: Change in temperature of the true (solid line) and the equilibrium (dashed 
line) transitions as function of t — t c ; = 30 GeV. 

Figure |9|: Change in temperature of the true (solid line) and the equilibrium (dashed 
line) transitions as function of t — t c ; ran = 60 GeV. 

Figure [l^: Change in temperature of the true (solid line) and the equilibrium (dashed 
line) transitions as function of t — t c ; ran = 80 GeV. 

Figure [IT]: Comparison of the temperature decrease due to supercooling and the tem- 
perature increase due to reheating as function of the Higgs boson mass. 



Figure [L2|: A two-soft-loop contribution to the effectiv potential. 
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Figure 1. 
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